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Abstract. We give an explicit formula for the joint density of the max and argmax 
of the Airy2 process minus a parabola. The argmax has a universal distribution which 
governs the rescaled endpoint for large time or temperature of directed polymers in 1 + 1 
dimensions. 



1. Introduction 

In geometric last passage percolation, one considers a family {w(i,j)} i j e %+ of inde- 
pendent geometric random variables with parameter q (i.e. F(w(i,j) = m) = q(l — q) m 
for m > 0) and lets II n be the collection of up-right paths of length n, that is, paths 
7r = (ttq, . . . , 7r n ) such that iTi — 7Tj_i G {(1, 0), (0, 1)}. The point-to-point last passage time 
is defined, for m,n G Z + , by 

m+n 

L(m,n) = max > w(-ir(i)), 

7ren m+n :(0,0)^(m,n) ^ 

where the notation in the subscript in the maximum means all up-right paths connecting 
the origin to (m, n). Next one defines the process t i— > H n (t) by linearly interpolating the 
values given by scaling L(n, y) through the relation 

L(n + y,n-y) = an + c 2 n 1/3 if n (c 3 n~ 2/3 y), 

where the constants q have explicit expressions which depend only on q and can be found 
in [Joh03]. The random variables 

T n = mf {t: supH n (s) = sup H n (s)} 

s<t seK 

then correspond to the location of the endpoint of the maximizing path with unconstrained 
endpoint. Johansson [Joh03] showed that 

H n (t) -+A 2 (t) -t 2 

in distribution, in the topology of uniform convergence on compact sets, where A2 is the 
Airy2 process, which is a universal limiting spatial fluctuation process in such models, and 
is defined through determinantal formulas for its finite-dimensional distributions (see the 
companion paper [CQR12] for a description). Together with known results for last passage 
percolation [BR01], Johansson's result (see also [CQR12]) implies that 

(1.1) p(sup{A 2 (t) - t 2 } <m)= F GOE (4 1/3 m), 

where -Fqoe is the Tracy- Widom largest eigenvalue distribution for the Gaussian Orthog- 
onal Ensemble (GOE) from random matrix theory [TW96] . 

Now let T denote the location at which the maximum is attained, 

T = argmax{„4 2 (£) — t 2 }. 

teR 

Together with the recent result of Corwin and Hammond [CHU] that the supremum of 
A 2 {t) — t 2 is attained at a unique point, Theorem 1.6 of [Joh03] shows 

1 
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Theorem 1. As n — > oo, T n — > T in distribution. 

In this article we complete the picture by providing an explicit formula for the distribu- 
tion of T ■ Let Ai denote the maximum of the Airy2 process minus a parabola 

M = max! A 2 (t) -t 2 }. 

teR 

Our main result is in fact an explicit formula (1.5) for the joint density f(t, m) of T and 
M. 

In the derivation of the formula, we will assume the result of Corwin and Hammond 
[CH11] that the maximum of A.z(t) — t 2 is obtained at a unique point. However, we point 
out that it is not necessary to do this. In fact, if one follows the argument without this 
assumption, one ends up with a formula for what is in principle a super-probability density, 
i.e. a non-negative function f(t,m) onlxl with J RxK dmdt f(t,m) > 1, and in fact one 
can see from the argument that 

/ dmdt f(t, m) = expected number of maxima of A%{t) — t 2 . 

JRxR 

Recall that from (1.1) that the distribution of Ai is given by a scaled version of -Fgoe- A 
non-trivial computation (see Section 3) on the resulting f(t, m) gives 



/oo 
dt/(t,m) = 4 1 / 3 ^ OE (4 1 /3 m) . 
-oo 



This shows that the resulting f(t,m) has total integral one, which can only be true if the 
maximum is unique almost surely. Thus we provide an independent proof of the uniqueness 
of the maximum of Ai (t) — t 2 . 

Now we state the formula. Let B m be the integral operator with kernel 

(1.2) B m (x,y) = Ai(x + y + m). 

Recall that Ferrari and Spohn [FS05] showed that -Fgoe can be expressed as the determi- 
nant 

(1.3) FgoeM = det(J - P B m P ), 

where P a denotes the projection onto the interval [a, oo) (the formula essentially goes back 
to [Sas05]). Here, and in everything that follows, the determinant means the Fredholm 
determinant in the Hilbert space L 2 (M). In particular, note that since -Fgoe^) > for 
all m G M., (1.3) implies that I — PoB m PQ is invertible. We will write 

Q m {x,y) = (I - PoB rn PQY 1 (x,y). 
Also, for t,m £M define the function 

(1.4) i(>t,m(x) = 2e xt [t Ai(x + m + t 2 ) + Ai'(x + m + t 2 )] 
and the kernel 

*t, m (x,y) = 2 1 / 3 ^, m (2 1 / 3 x)^_ t , m (2 1 / 3 y). 

Finally, let 

rOQ roc 

7 (t,m) = 2 1 / 3 / dx dy^_ tA „ 1/3m (2 1 / 3 x) gm (x,y)4> tA - 1/Sm (2 1 / 3 y). 
Jo Jo 

Theorem 2. The joint density f(t,m) ofT and Ai is given by 

f(t, m) = j(t, 4 1 / 3 m)F GO E(4 1/3 m) 
1 = det(J - PoB^, 3m Po + PoVt,mPo) ~ F GO E(4 1/3 m). 
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Integrating over m one obtains a formula for the probability density / en d(^) of T ■ Un- 
fortunately, it does not appear that the resulting integral can be calculated explicitly, so 
the best formula one has is 

/oo 
dmf(t,m). 
-oo 

One can readily check nevertheless that f en d(t) is symmetric in t. In [QR12] it is shown 
that the tails decay like e~ ct . Figure 1 shows a contour plot of the joint density of Ai and 
T, while Figure 2 shows a plot of the marginal T density. The numerical computations 
of Fredholm determinants used to produce these plots are based on the numerical scheme 
and Matlab toolbox developed by F. Bornemann in [BorlOa; BorlOb]. 

Although one only has the rigorous result in the case of geometric (or exponential) last 
passage percolation, the key point is that the polymer end-point density f en d(t) 1S expected 
to be universal for directed polymers in random environment in 1 + 1 dimensions, and 
even more broadly in the KPZ universality class, for example in particle models such as 
asymmetric attractive interacting particle systems (e.g. the asymmetric exclusion process), 
where second class particles play the role of polymer paths. And the analogous picture is 
expected to hold, as we now describe. 

In the directed polymer models we consider a family {w(i, j)}«ez+ j e z of independent 
identically distributed random variables and the probability measure (polymer measure) 
P™ a on the set Tl n of one-dimensional nearest-neighbor random walks of length n starting 
at given by 

where f3 > is the inverse temperature. The analogue of T n in this context is 7r(n), the 
random position of the endpoint. The last passage percolation case corresponds to f3 = oo. 
The infinite temperature case (3 = is nothing but a free random walk. For (3 < oo the 
endpoint is random even given the random environment {w(i, j)}jez+jez- Still one expects 
in great generality, and for any /3 > 0, to have 

cn~ 2/3 ir(n) / end 

for an appropriate c. The conjecture is that this holds whenever < oo, and fails 

otherwise due to the appearance of special large values of w which attract the polymer. 
However, few results are available at finite temperature. The first model for which any 
results were obtained (for the free energy) is the continuum random polymer (see below). 
There are now two other models, the semi-discrete model of O'Connell-Yor [OY01], and 
the log-Gamma polymer [Sepl2; COSZ11], for which results about asymptotic fluctuations 
of the free energy are becoming available. 

In the context of the continuum random polymer, we have continuous paths x(s), < 
s < t, starting at at time 0, with quenched random energy 

nx(-)) = ft{\x( S )\i-Z(s,x(s))}ds, 

where £ is Gaussian space-time white noise, that is, (£(i, x),£(s, y)) = 5(t — s)5(y — x). 
Through a mollification procedure [AKQ12] one can construct a probability measure Pf 
on the space of continuous paths corresponding to the formal weights e - ^. It has finite 
dimensional distributions P^{x{t\) G dxi, . . . ,x{t n ) G dx n ,x{t) G dx), < t\ < ■ ■ ■ < t n < 
t, given by 

Z(0, 0, t\, x\) ■ ■ ■ Z (t n -i, x n _i, t n , x n )Z(t n , Xfi ,t,x) 

7^- ; dx\ • • • dx n ax 

JdyZ(0,0,t,y) 

where Z(s, y, t, x) is the solution of the stochastic heat equation with multiplicative noise 

d t Z = p~ l d 2 x Z + PiZ 
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Figure 1. Contour plot of the joint density of M. and T ■ 



on (s,t] with initial data Z(s,y,s,x) = 5(x — y). The temperature can be related to time 
as t ~ /3 4 , so through a time rescaling we can set (3 = 1 without loss of generality. 

In this setting the endpoint distribution is 

Z(0,0,t,x) 



P?(x(t) G dx) 



fdyZ(0,0,t,y) 



dx. 



Writing 
(1.7) 



Z(0,0,t,x) 



+{UY^A t {{U)- 2 ^x)+ 



21 



oo, the crossover process 



the key prediction (see Conjecture 1.5 of [ACQ11]) is that, as t 
At converges to the Airy2 process, 

A t {x) -^A 2 {x). 

This is proved in the sense of one dimensional marginals in [ACQ11; SS10], and a non- 
rigorous computation for multidimensional distributions was made in [PS 11]. Calling x = 
(4i) -2 / 3 x we can rewrite the exponent in (1.7) as (4t) 1 ^ 3 {At(x) — x 2 } + X, from which we 
conclude that the endpoint of the polymer at time t has approximately the distribution 
(4£) 2 ' 3 7" for large t. The partition functions of discrete directed polymer models satisfy 
discrete versions of the stochastic heat equation, and analogous results are expected to 
hold in that setting as well. 

The problem has attracted interested in the physics literature for quite some time (see for 
example [MP92; HHZ95]). Recently there has been a resurgence of interest. In particular, 
an alternate way to obtain the Airy2 process is as a limit in large N of the top path in 
a system of N non-intersecting random walks, or Brownian motions, the so called vicious 
walkers [Fis84]. Schehr, Majumdar, Comtet, and Randon-Furling [SMCRF08], Feierl [Fei09] 
and Rambeau and Schehr [RS10; RS11] obtain various expressions for the joint distributions 
of A4 and T in such a system at finite N. Forrester, Majumdar, and Schehr [FMS11] obtain 
the -Fgoe distribution from large N asymptotics non-rigorously, and furthermore make 
connections between these problems and Yang-Mills theory. Unfortunately, the formulas 
obtained for T at finite N have not been amenable to asymptotic analysis. 
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Figure 2. Plot of the density of T compared with a Gaussian density with 
the same variance 0.2409 (dashed line). The excess kurtosis E(7~ 4 ) /E(T^) 2 — 
3 is -0.2374. 
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2. Derivation of the formula 

Let (Ml,Tl) denote the maximum and the location of the maximum of .42 (t) — t 2 
restricted to t E [— L,L], and let /x be the joint density of (M.l,Tl)- We first note that, 
by results of I. Corwin and A. Hammond [CH11], the joint density f(m, t) of A4,T is well 
approximated by /i(m, t), 

(2.1) f(t,m) = lim f L (t,m). 

L— >oo 

By definition, 

f L (t,m) = lim lim^-P(A4 L E [m,m + e], T L G [t,t + 5]), 
<5^0 £-s>0 SO 

provided that the limit exists. The main contribution in the above expression comes from 
paths entering the space-time box [t, t + S] x [m, m + e] and staying below the level m 
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outside the time interval [t, t + 6]. More precisely, if we denote by D e s and D £i s the sets 
D eS = {a 2 (s) - s 2 < m, s G [t,t + 5] c , A 2 (s) - s 2 <m + e, s€[t,t + S\, 

A 2 (s) — s 2 G [m, m + e] for some s G [t,t + <5] |, 

and 

D £ ^s = ^A 2 {s) — s 2 < m + e, s G [— L, L], .42 (s) — s 2 G [m, m + e] for some s G [t, t + 6]\, 
then 

D e S C {X L G [m,m + e], 71 G [t, t + 5]} C 7J £)(5 . 

Letting f(t,m) = lim^o nm e->o ^^(iZe^) and defining f(t,m) analogously (with D £ ^ 
instead of De s) we deduce that f(t,m) < f(t,m) < f(t,m). In what follows we will 
compute f(t, m). It will be clear from the argument that for f(t, m) we get the same limit, 
so we will only compute f(t,m). The conclusion is that 

/ L (t,m) = hmhm^-P(/2 £i5 ). 

S— >o e— >o so 

We rewrite this last equation as 
(2.2) 

f L (t,m) = lim lim-JPGAaOO < h E , 5 {s), s G [-L, L]) - F{A 2 (s) < h , S (s), s G [-L,L])], 

5-5-0 e->0 £0 1 J 

where 

h s ,s(s) = s 2 + m + el se[t)t+S y 

Our method is based on precise computation of the two probabilities. We recall the formula 
in Theorem 2 of [CQR12] for the probability that A 2 (t) < g(t) on a finite interval. Introduce 
the operator ©j^ r ] which acts on L 2 (M) as follows: ®fg r i/(") = u ( r i')i where u(r, •) is the 
solution at time r of the boundary value problem 

d t u + Hu = for x < g(t), t G (£, r) 
u(£,x) = f(x)l x<g{e) 
u(t, x) = for x > g(t) 

for the Airy Hamiltonian, 

H = -d 2 x + x. 

In [CQR12] it is shown that this operator describes the height statistics of the Airy2 process, 
(2.3) F(A 2 (t) < g(t) for t G [£, r}) = det(7 - K M + e~ iH K A M ,e rH K M ) , 



where we have used the cyclic property of determinants as in (1.7) in [CQR12]. We use 
(2.3) to rewrite (2.2) as 



1 



f L (t, m) = lim lim - det ( J - K Ai + e LH K Ai Q^[ L] e LH K A 



S-s-0£->o eS - 



The limit in e becomes a derivative 

1 



det (/ - K A [ + e LH K Ai Q h ^ s L L] e LH K Ai 



f L (t,m) = lim - dp det [I - K Ai + e hH K Ai Q'^ L] e hH K Ai 



/3=0 
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which in turn gives a trace, 



2.4) f L (t, m) = detf I - K Ai + e LH K Ai G'£l L] e LH K A 



i. h o,s LH 



lim — tr 

8^0 



(I - K Ai + e LH K Ai e'(y LM e^K Ai y L e^K A 



ho,S LH jy \-lLH 



hp, 6 
-L,L] 



13=0 



e LH K M 



(see Lemma A. 2 and Remark A. 3). Note that hog = 9m, where g m is the parabolic barrier 



g m {s) = s +m, 

so in particular the determinant and the first factor inside the trace do not depend on 5. 
Prom (1.2) and Theorem 1.3 from [CQR12] we have 



(2.5) 



ton (I - K Ai + e LH K M Q^ s L L] e LH K Ai j =1- AP R.P Q A* 



in trace norm, where P a = I — P a denotes the projection onto the interval (— oo, a], 

B}(X, A) = 2- 1 / 3 Ai(2- 1 / 3 (2m — A — A), 
and the Airy transform, A, acts on / £ L 2 (JR) as 



Af(x) 



dz Ai(x — z)f(z). 



In particular, (1.8) in [CQR12] implies that 



(2.6) 



lim det 

L— >oo 



I-K Ai + e LH K A ^ 6 LM e hM K Ai \ = F GOE (4 1/a m 



x^o,«5 LH 



The next step is to compute dpQ,^'£ L , |^=o- Recalling that ho^(s) = g m (s) = s 2 + 
and also h e< s(s) = g m+e (s) for s £ [i, t + <5] we have, by the semigroup property, 



m 



[—L,L] [—L,L] ~ [-L,t] 



C\9m+s p\9m 



ClSra 



We now use Theorem 3 of [CQR12] and a minor variation of (1.4) in [CQR12] to obtain 

that ©njyi nas explicit integral kernel 

(2.7) 

e £x-ry+(r 3 -e 3 )/3 



4(r-£) 



4(r-£) 



c<m+£ 2 *-y<.m+r 2 • 



For convenience we introduce the kernels #i(x, z) = e t2: 0r^ ^(sc, z)l x < m +L2 and $2(z, y) = 

e _<z 0^^ L j (z, y)~Ly<m+L 2 > where ©r/j is defined as in (2.7) but with the indicator functions 
replaced by 1. Let 



(2.8) A^(x,y) 



m+t 2 



dz 



m+(t+5) 2 



dz'&ilx, z) 



e -(z-t 2 +z~(t+6) 2 -2m) 2 /(45) _ e -(z-t 2 +S-(t+S) 2 ~2m-2e) 2 /(AS) 



which corresponds to l] ~®[~l l] without shifting m by e in the indicator functions 
in (2.7) for the first operator in this difference. We will show in Lemma A. 4 that 



(2.9) 



lim - 

£->0 e 



V{-L,L] [—L,L] 



A 



e,<5 







in Hilbert-Schmidt norm. On the other hand, performing in (2.8) first the change of 
variables zi-yz + m + f 2 , z, z + m + (t + 5) 2 , then a scaling of z and z by Vo 5 and then 
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the change of variables — u = z + z, —v = z — z,we get 

e [(t+5) 3 -t 3 ]/3 poo f-u 

a l ( x > v) = - 



4v^ 



du I dvtf 1 (x,-V5(u + v)/2 + m + t 2 ) 

2 (Vtf(v - «)/2 + m + (t + 5) 2 , y), 



e -n 2 /4 _ e -(u+2e/VS) 2 /i 



From this form and (2.9) it is straightforward to see that 



(2.10) lim - 

£->o e 



[—L,L] ~ [—L,L] 



(x,y) = lim -A £ jf(x,y) 



4Jn 



OO CM 



ciu / dvue~ u /4 $ 1 (x,-V5(u+v)/2+m+t 2 )$ 2 {V8(v-u)/2+m+(t+5) 2 ,y). 



The limit holds in Hilbert-Schmidt norm, as will be shown in Lemma A. 4. Now we take 
the limit in 5 and obtain 



(2.11) 



lim - 

<5-s-o 5 



/3=0 



(x,y) = d w -&i(x,w)\ w=m+t 2 d w i} 2 (w,y)\ w=m+t 2, 



again in Hilbert-Schmidt norm, which will be checked in Lemma A. 4. Referring back to 
(2.4) we have now shown that 

(2.12) " ' 



lim - x e LH K Al 
5^0 



where Vl/ l has kernel 

with 

(2.13) 



0=0 

V L {x,y) = ^>l{x)^ 2 L (y) 



e LH K Ai = * L , 



®l(x) = d w (e LH K Ai Qp Lt] M t (x,w)) 
*l(y) = d w (M- t @f™ L] e LH K Ai (w,y)) 



w=m+t 2 

w=rn+t 2 
tx . 



and Mt is the multiplication operator given by Mtf(x) = e f(x). 

Putting (2.1), (2.4), (2.6) and (2.12) together and using Lemma A. 1(a) we have 



(2.14) f(t, m) = lim tr (/ - K Ai + e LH K A] Q 9 r T n e Ltl K M )- l V L F GOE (4 1/a m). 



lV3, 



L— >oo 



We now have to compute the limit of the trace. We begin by using (2.7) to compute 



ip(z) :=d w [Qf_ LA M t {z, W )) 



c -Lz+L 3 /3+t 3 /3 

= e -—, -^{z -m-L 2 ) e -(*—i a ) a /4(£+t). 

Note how the derivative of the two terms inside the bracket in (2.7) evaluated at w = m + t 2 
are equal. From (2.13) we get 

V\{x) = e LH K Ai P m+L2 ip{x) = e LH K Ai ip(x) - e LH K Ai P m+L 2 tp(x), 
In Appendix A we will show that 

(2.15) ljm \\e LH K Ai P m+L 2 <p\\ L2{R) = 0. 

Now we compute e LH K A np. We write it as 



(2.16) 



e LH K Ai ip(x)= I dX e Alj Ai(x — A) / dz Ai(z - A) ip(z) 



XL 



To compute the z integral, which we denote by /(A), we use the contour integral represen- 
tation of the Airy function given by 



(2.17) 



Ai(x) = — I due u3/3 - ux , 
2-ki Jv 
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with r = {c + is : s E M} and c any positive real number, to write 



„-L(z+m+L 2 )+L 3 /3+t 3 /3 
/(A ) = ^ I du I dz e ^/Z-u( Z+ m + L*-X) -* KL -, , 



2^{L + t)3/2 



1 

2vri 

where we have shifted the variable zbym + L 2 . Note that the integral in z is of the form 
dz e~ aiz ~ a 2 z ~ a s z ^ w hi c h corresponds to computing the mean of a certain Gaussian 
random variable. Performing the integration we get 



/(A) 



2vri 



J du e " 3 /3+(i+<)« 2 + (L 2 +2Lt-m+A)u-Lm+L 3 /3+L 2 t+t 3 /3^ _|_ 



U . 



Introducing the change of variables u = v — L — t we get 

/(A) = - A e ^ 3 -(L+t)A /" due „3 /3 _ (m+t 2_ A)v ^ _ 
27Ti 7 r / 

where T' corresponds to a shift of T along the real axis. Using (2.17) we deduce that 
/(A) = 2 e mt+t 3 -{L+t)\ r Ai /( m _|_ _ A) + f Ai(m + t 2 - A) . 

Therefore 

/•o r 
e LH K M <p(x) = 2 dX e * 3 +( m - A )* Ai(x - A) Ai'(m + t 2 - A) + t Ai(m + i 2 - A) 

We will rewrite this identity as 

e LH K Ai cp(x) = AP $ t>m (x), 



where 



Ai'(m + i 2 - x) + t Ai(m + i 2 - x) 



Remarkably, the result does not depend on L. Note that APQtpt,m G -^ 2 (^); which can be 
checked using the Plancherel formula for the Airy transform 



(2.18) 



{Aff 



f 



and the fact that |Ai(tt)| V \Ai'(u)\ < Ce~h s/2 for some C > and all u > (see (10.4.59- 
60) in [AS64]). 

Now we look at ^\{y). By the time symmetry and time homogeneity of the heat kernel 
it is clear that d w (M- t @?™ L] (w, -))(y)\ w= . i9 can be obtained from the above calculation 
by starting at y and running backwards in time from L to t. Observe that the length 
of this time interval is L — t, whereas the one in the above calculation had length L + t. 
Moreover, here we are multiplying the boundary value operator by M-t, whereas before 
we multiplied by Mf. It is not difficult then to see that the answer for the second factor 
should be the same as for the first one, only with x replaced by y and t by —t. From this, 
(2.13) and (2.15) we get that 

V L (x,y) ► V(x,y) := AP $ t!rn (x)APoi)-t,m(y) 

L— >oo 

in Hilbert-Schmidt sense, and thus from (2.5) and Lemma A. 1(b) we have that 



(/ - K Ai + e LH K Ai Q^l L] e LH K Ai )- L ^L ► (I - AP R'P A 
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in trace norm (the product converges in trace norm thanks to Lemma 3.1 of [CQR12]). 
Therefore by Lemma A. 1(a), 



lim tr 

L— >oo 



(/ - K M + e LH K M Q^l L] e LH K M )- l ^> L 
tr 



-L,L 
?1 f3 a*\— 1 



(/ - AP R L P A*)- L V = ((I - AP R L P A*)- x APoflH,m,AI\ii/)-t, m . 

where (•,•)% denotes inner product in the Hilbert space V. (with U = L 2 (R) if the subscript 
is omitted). 

It only remains to simplify the expression. We use the reflection operator crf{x) = f{—x). 
Because (Aa) -1 = Aa, a 2 = I and A* = aAa, we have 

(I - AP R 1 PoA*y 1 AP i, t , m ,APo^-t, m ) = (a<t{I - aP R l P Q a)- l aP^t >m ,AP^ t>r 
Since (Aa)* = Aa and AaA = a, this last term can be rewritten as 

(I - a/WiVrVPoVVn, (Aa)*AP 4>-t, m ) = ((I - aP Q R}P Q a)- x aP^ t , m ,aP^ t , m 

(I - P aR l aPo)~ 1 ai!t,m,cri'-t,m) , 

/ L 2 ([0,oo)) 

where in the second equality we used the trivial fact that Poa = aPo and aPo = Poa. 
Observing that aipt,m(x) = e* +mt ipt,m(x), where ipt,m was defined in (1.4), we deduce that 

((I - AP G R l P Q A*)- l AP^ tym ,AP Q i,_ tym ) oi = ((I - P Q aE}aP Q )-^ t ,mA-t, m ) on „ 

\ / L 2 (R) \ I L 2 ([0,oo)) 

Now we use the scaling operator Sf(x) = /(2 1 / 3 x). One can check easily that S^ 1 = 
2 1 / 3 S* and that P commutes with S and S' 1 . Since aR l a(x,y) = 2 -1 / 3 Ai(2" 1 / 3 (a; + y) + 
4 1 / 3 m), we also have 

SanaS 1 = i3 4 i/3 m , 
where this last kernel was defined in (1.2). Thus writing m = 2~ 1 / 3 m we get 

({I - PoO--R 1 (TPo) _1 ^,m,V'-i,m^ 2([0oo)) = ((I ~ S~ 1 PoB 2 fhPoSy 1 ^t,m,ij-t,m) L 2( [0i0c y ) 

= (S^il ~ PoB 2 rhPo)' 1 SA,m,^-t,m) L 2 {[0tOD)) 
= 2 1 / 3 ((J - P B 2 ^P )- 1 S^ m , ^-t,m) i 2 ([0)Oo)) , 

which is equal to 2 1 / 3 7(t, 4 1 / 3 m). This gives our first formula for f(t,m) in (1.5). Now 
observe that 7(t,4 1 / 3 m) equals the trace of the operator (/ — PQB i i/3 m Po)~ 1 Po^ t ,mPo 
and that ty t m is a rank one operator. The second equality in (1.5) now follows that 
from the general fact that for two operators A and B such that B is rank one, one has 
det(J - A + B) = det(J - A) [l + tr ((/ - A)~ 1 B)} . 

3. M MARGINAL AND UNIQUENESS OF THE MAXIMIZER 

As we mentioned in the introduction, Corwin and Hammond [CH11] showed that the 
maximum of A2(t) — t 2 is attained at a unique point iGR, providing a proof of a conjecture 
by K. Johansson (Conjecture 1.5 in [Joh03]). We used their result in Section 2 to write 
formulas for f(t, m) in terms of certain events concerning the Airy2 process. 

Alternatively, one can turn the reasoning around and use our formula to give a different 
proof of Johansson's conjecture. If we do not assume the uniqueness of the maximizer, 
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then the derivation in Section 2 leads to a density f(t,m) for the event that there is a 
maximizer at t (and height m). Therefore the uniqueness of the maximizer is equivalent to 

/ dt dm f(t, m) = 1. 
Jm 2 

This, in turn, is a direct consequence of the following 
Proposition 3.1. For any m € M, 

/OO J 
dtf(t,m) = — F GOE (4 1/3 m). 

Proof. From the formula (1.5) for f(t,m) we see that we need to compute 

POO 

W m (x,y)= / dtlj}-t,m{x)i>t,m{y), 



where x = 2 1//3 x and y = 2 1 / 3 y. Let r a = {a + is : s € M}. Then fixing a > and using 
(2.17) we have 

4 



(4m) 2 



dudv / dt (u — t)(v+t)e 



u 3 /3+v 3 /3-u(x+m+t 2 )-v(y+m+t 2 )+t(x-y) 



The t integral is just a Gaussian integral and gives 



(4vri) 2 



r 2a xr Q 



dudv (it + t;)- 5/2 ^, s (n, uje* 8 '*^, 



where 

Px, y (u, v) = Av?v + 4uv 3 + 8u 2 u 2 - 2 (it + i>) + 2(u 2 - w 2 )(x - y) - (x - yf 
and 

qx,y(u,v) = - 



-(it 4 + 1> 4 + it 3 i; + uv 3 ) — m(u + v) 2 — u 2 x — v 2 y + \{x — y) 2 — uv(x + y) 



u + v 

Introducing the change of variables z = u + v, w = u — v, we get 

R 1 



^m{x,y) 



(4vri)2 2 J TAa J Ta 



dz / dwz- 5 / 2 p £ ^(z,ii;)e^^ (2 ' ,i ' ) , 



where 



and 



Px,y(z, w) = —w z + 2wz(x — y) — (x — y) — 2z + z 



w 2 z 2 — 2wz(x — y) + (x — y) 2 — 2{x + y + 2m)z 2 + \z A 

Iz 



qx,y(Z,W) 



Changing variables w i— )■ iio, the it; integral is another Gaussian integral and we get 



^ m (x,y) 



1 
47ri 



dzze z 3 /12-z(x+y+2m)/2 



4 2/3 

47ri 



dz z e z 3 /z-z- 1/A z{x+y+2rh) 



r 4-!/3 3a 



= -2 1 / 3 Ai , (x + y + 4 1 / 3 m)) = -2^d m B Al/Zm (x, y), 
where we have used (2.17). Using this in the definition of j(t, m) we deduce that 



dt-y(t,m) = -2 l 'Hr (/ - P Q B m P Q )- 1 * r 
Consequently we get from (1.5) and (1.3) that 



tr 



(I — Po-Bm-Po) 1 d m B Tl 



dtf(t,m) = -tr[(J - PoBmPo)' 1 d m B m \ det(l - P PmP ) = ^ det(/ - P P m P ), 
where the last inequality follows from Lemma A. 2. The result now follows from (1.3). □ 
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Appendix A. Technical estimates 

Section 3 of [CQR12] contains a short review of some general facts about trace class and 
Hilbert-Schmidt operators and Fredholm determinants. In Section 2 of the present article 
we used some additional facts, which we state next. Here H will denote a separable Hilbert 
space and Bi(H) will denote the space of trace class operators in H, which is endowed with 
the trace norm (see Section 3 of [CQR12] for a short discussion or [Sim05] for a complete 
treatment). 

Lemma A.l. Assume {A(v)} >Q is a family of operators converging as v —> oo in £>i("H) 
to some operator A G B\(7-l). Then: 

(a) ti(A(v)) >ti(A). 

V — ^oo 

(b) If I — A(v) is invertible for all large enough v and I — A is also invertible, then 

(I-A(v))- 1 >(I-A)- 1 inBiCH). 

v— >oo 

This result comes from Theorem 3.1 and Corollary 5.2 in [Sim05]. Using (5.1) from 
[Sim05] one can also easily show the following (see also the corollary just cited): 

Lemma A. 2. Assume \A(f3)}p^ ^ is a family of operators in B\{H) such that there is 
an operator dpA{<S) satisfying 

l[A(J3)-A(0)]j-+d f) A(fi) inBxCH). 

Then the map f3 i — > dei(l + A{j3)\ is differentiate at and 

dp det(J + A(p)) = tr[(J + ^(O))- 1 ^^)] det(J + A(0)). 

Remark A. 3. Note that the last two lemmas assume convergence in trace norm as the 
hypothesis. Throughout Section 2 (see (2.4), (2.10) and (2.11)) we used these results for 
operators of the form e LH K A ^ rj e LH K A i, where r\ is some parameter and we know that 
^ converges in Hilbert-Schmidt norm to some limit <£. As we will see in Lemma A. 4, 
the convergence is in fact a bit stronger, and using this we can justify the application of 
the lemmas in Section 2. To see why, note that if we let (p(x) = 1 + x 2 and define the 
multiplication operator Mf(x) = ip(x)f(x) then by Lemma 3.1 of [CQR12] we have 

(A.l) \\e LH K M {% - ^e^Ailli < ||e LH ^ A i||op||(^ - $)M|| 2 ||M- 1 e L ^ A i||2. 
By (2.16) we have, for / G L 2 (R), 

\\e LH/2 K Ai f\\l = I dxdydyf d\ dX e ^ L / 2 Ai(x - A) M{y - X)f(y) 

JR 3 J(-oo,0] 2 

•Ai(x-A)Ai(y-A)/(y) 
dy dy [ dX dX e ^ L ' 2 Ai(y - X)f{y) Ai(y - X)f(y)5 x= ~ x 

J(-oo,0] 2 

dAe AL ,4/(A) 2 . 

Using (2.18) we deduce that ||A|| p = ||^4*|| p = 1, and then 
(A.2) \\e LH / 2 K Ai \\ op < 1. 

The third norm in (A.l) is also finite, thanks to (3.3) in [CQR12], and we are going to 
prove below the convergence ||(4>ti — < 1 ) )M||2 — > in each relevant case. 



The next result provides the missing estimates in the proof of (2.14). 



ENDPOINT DISTRIBUTION OF DIRECTED POLYMERS 



13 



Lemma A. 4. For each fixed 5,L > 0, the convergences in (2.9), (2.10) and (2.11) hold 
in Hilbert- Schmidt norm. Moreover, if we let <p(x) = 1 + x 2 and define the multiplication 
operator Mf(x) = ip(x)f(x), then the three convergences above still hold if we multiply 
each operator on the right by M. 

Proof. The second equality in (2.10) one follows from the dominated convergence theorem 
and the estimate 



£ 



e « 2 /4 _ e -(u+2e/V5) 2 /A 



ue 



<cj=(l + u 2 )e- u2 /\ 



where C > can be taken uniform in u > for small enough e. Using this bound and the 
particular form of i?i and "Q 2 we can see that 



du I dv 



e -« 2 /4 _ e ~(u-2e/VS) 2 /A 



ue 



■ &i(x, -Vd(u + v)/2 + c + t 2 ) # 2 (y/6(v - u)/2 + m + (t + 5) 2 , y) 



< C 



M\x\+\v\)- 



for some C > 0. Integrating the square of the left side with respect to x and y over 
(— 00, m + L 2 ] 2 , we can deduce again by the dominated convergence theorem that e~ 1 A e £ S 
converges in Hilbert-Schmidt norm. This, together with (2.9), proves (2.10). 

Next we observe that 

#i(x, -VS(u + v)/2 + m + t 2 )$ 2 (V5y(v -u)/2 + m + (t + 5) 2 ,y) 

+ j(u + v)(v - u)d w $i(x,w)\ w=m+t 2 d w $ 2 (w,y)\ w=m+t 2 < 5 i/2 e(u, v, x, y), 

where e involves products of first and second derivatives of i?i and v x 2 . By the same argu- 
ment we explained above, J du J dv \e(u,v, x,y)\ can be easily seen to be in L 2 {{— 00, m + 
L 2 ] 2 ) as a function of x and y. Thus, by the dominated convergence theorem, 



lim - 

<5->o 



1 



[-L,L 
00 

du 



13=0 







U 2 1 

dvu{u + v)(v - u) e~ u2/i -d w v'i(x,w)\ w=m+t 2 d w -d 2 (w,y)\ 

w=m-\-t 2 

-11 ^ 



in L 2 ((— 00, m + L 2 ] 2 ). The integral in u and v can be computed, and gives the answer 
— 16-v/7r, so we deduce that 



1 

6 L 



hp.s 



(x,y) - — > d w $i(x,w)\ w=m+t 2 d w $ 2 (w,y)\ w=m+t 2 
p=0 0—^0 



in Hilbert-Schmidt norm. This proves (2.11). 

We are left with proving (2.9). Let E e = (©j^ £] — l]) ■ ^° simplify notation 

we assume m = t = 0, for the general case the proof is exactly the same. Prom (2.7) and 
(2.8) we have 



E s {x,y) 



53/3 / dzdzMx,z) 



D 



e -(z-z+8 2 f/(A5) _ e -(z+z-S 2 -2s) 2 /(A5) 



where D = ((— 00, e] x (— 00, e + 5 2 ]) \ ((— 00, 0] x (— 00, 5 2 ]). We split D into the union 
of three disjoint regions of pairs (z, z): D\ = {0 < z < e, 5 2 < z < 5 2 + e}, D 2 = {0 < 
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z < e, z < 0} and D% = {z < 0, 5 2 < z < 5 2 + e}. Similarly we split E e as the sum of the 
integrals El over each region. On the first region we have 

- \E l e (x,y)\ < -^=e s3/3 - [ dzdzMx,z)#2(z,y) >® 



'4x5 e J Dl e^o 

thanks to the particular form of $i and $2 and the fact that D\ has area e 2 . For the second 
region we have 



= 0, 

2=0 



lim -E z Ax, y) = ^= / dz&i(x, z) e -i--^ 2 ) 2 /m _ e -(z+z-8^/(46) y) 

£^0 £ V47T<5 J-oo L J 

while the third region can be dealt with analogously. We deduce by the triangle inequality 
that e~ 1 \E £ (x,y)\ — > as e — > 0. To upgrade the convergence to Hilbert-Schmidt norm we 
may use the dominated convergence theorem and similar estimates as for (b) and (c), we 
omit the details. This finishes the proof of (a). 

Finally, it is straightforward to check in each case that the convergences still hold if we 
multiply each kernel by the polynomial 1 + y 2 . □ 

Proof of (2.15). By Lemma 3.1 and (A. 2) we have 

\\e LH K Ai P m+L 2(p\\ 2 < He^KAillopH-Pm+L^lh < \\P m +L 2( Ph- 
This last norm can be easily computed: 

„ 9 1 2,3, 2,3 f°° , 9n9 ( Z -m-L 2 ) 2 2L 

il^m+^^ill^ ^ ,t , ^ 5 / dz(z-m + L 2 ) 2 e ^+*) 



167r(L + t) 3 Jm+L 2 



! r e-! i3 +I i3 f dz(z- m ) 2 e-^- 2Lz . 



16vr(L + t) 3 

Let Fl(z) denote the argument in the last exponential. Fl is minimized at z* = m — 
2L{L + 1), which is less than m for large L, and is strictly increasing in [z*,oo). Thus Fl 
attains its minimum inside the interval [m, oo) at z = m, where its value is —2mL. An 
application of Laplace's method (Lemma 5.1 of [CQR12]) then shows that 

\\Pm+LM\l < Ce~ L3 / c , 
for some C > 0, which finishes the proof. □ 
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